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The  analysis  of  the  electrochemical  behaviour  of  polymer  electrolyte  fuel  cells  (PEFC)  both  in  time  and 
frequency  domain  requires  appropriate  impedance  models.  Simple  impedance  models  with  lumped 
parameters  as  resistances  and  capacitances  or  Warburg  impedances  do  have  limitations:  often  the  validity 
is  limited  to  a  certain  frequency  range,  effects  at  very  low  or  very  high  frequencies  can  not  be  described 
properly. 

However,  these  models  have  their  usefulness  for  engineering  applications,  e.g.  to  distinguish  the  major 
loss  terms,  to  estimate  the  membrane  resistance,  and  to  observe  the  changes  of  internal  losses  of  fuel  cells 
over  time  without  the  need  for  additional  sensors. 

The  work  discusses  different  impedance  configurations  and  their  applicability  to  impedance  spectra  of 
a  fuel  cell  stack.  Impedance  spectra  at  points  along  the  DC  polarization  curve,  as  well  as  spectra  at  various 
operating  conditions  are  analysed  and  identified  by  a  complex  nonlinear  least  squares  method.  Finally, 
the  connection  of  the  impedance  data  with  and  the  assignment  of  the  parameters  to  physical  phenomena 
are  discussed.  The  examination  shows  that  simple  impedance  models  are  well  qualified  to  describe  the 
electrochemical  behaviour  over  a  wide  frequency  range  at  all  operating  conditions. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  static  IV-characteristics  of  a  polymer  electrolyte  fuel  cell 
shows  a  nonlinear  behaviour.  The  major  losses  occurring  in  the 
electrochemical  reaction  are  due  to  the  charge  transfer  through 
electrode-electrolyte  interfaces,  the  transport  of  gases  through  the 
porous  layer  and  the  conduction  of  protons  through  the  polymer 
electrolyte  membrane.  The  losses  and  their  dependency  on  the  cur¬ 
rent  density,  as  depicted  schematically  in  Fig.  1,  are  theoretically 
described  and  modelled  in  terms  of  overvoltages  [1,2]. 

Nonetheless,  a  distinction  of  the  single  processes  and  their 
contribution  to  the  overall  losses  is  not  possible  by  static  measure¬ 
ments.  IV-characteristics  do  not  provide  an  insight  in  the  fuel  cell.  To 
describe  the  dynamic  electrochemical  behaviour  of  a  fuel  cell  and 
to  distinguish  single  loss  factors,  further  information  and  dynamic 
measurements  are  needed. 

For  the  characterisation  of  the  current  state  of  a  fuel  cell,  measur¬ 
able  states  as  pressures,  temperature,  and  mass  flows  are  available 
which  are  measured  at  the  inlet  or  outlet  of  the  fuel  cell.  Further¬ 
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more,  the  load  current  and  the  resulting  fuel  cell  voltage  are  known. 
Internal  loss  terms  which  are  crucial  for  the  resulting  voltage  and 
therewith  for  the  efficiency  of  the  fuel  cell  cannot  be  measured 
directly.  Hence,  for  the  operation  and  development  of  fuel  cells 
methods  are  required  which  provide  information  on  internal  losses 
and  thereby  enable  an  analysis  of  the  processes  taking  place  inside 
the  fuel  cell.  For  the  operation  of  fuel  cells  information  on  the  gas 
supply,  on  the  humidification  of  the  proton  conducting  membrane, 
and  on  transport  hindrances  are  needed.  For  the  production  of  fuel 
cells  occurring  ohmic  resistance  of  membranes  and  electron  con¬ 
ducting  components  need  to  be  measured  and  validated.  For  the 
development  of  fuel  cells  processes  leading  to  degradation  need  to 
be  analysed. 

An  approach  which  gives  an  insight  in  the  fuel  cell  stack  is 
the  analysis  and  identification  of  the  fuel  cell  impedance  [3]. 
The  impedance  ZFC(j&>)  can  be  regarded  as  the  complex,  nonlin¬ 
ear,  frequency-dependent  internal  resistance  of  the  fuel  cell.  The 
impedance  of  a  fuel  cell  incorporates  detailed  information  on  the 
intrinsic  loss  factors,  on  the  conductivity  of  the  membrane,  on  the 
electrode  processes,  and  on  kinetic  losses  [4].  The  measurement 
of  the  impedance  spectrum  at  an  operating  point  and  the  iden¬ 
tification  of  an  appropriate  impedance  model  enable  the  desired 
distinction  of  loss  terms,  and  a  description  of  the  dynamic  electro¬ 
chemical  behaviour  of  a  fuel  cell. 
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Fig.  1.  IV-characteristics  of  a  polymer  electrolyte  fuel  cell  with  charge  transfer  over¬ 
voltage  ijct,  ohmic  overvoltage  and  transport  overvoltage  ^trans- 

The  most  widely  used  method  for  analysing  impedances  is  the 
electrochemical  impedance  spectroscopy  (EIS).  The  approach  of 
EIS  is  to  measure  impedances  by  applying  successively  single¬ 
frequency  currents  to  the  system  and  measuring  the  real  and 
imaginary  parts  of  the  resulting  voltage  at  that  single  frequency. 
Measuring  the  frequency  response  yields  an  impedance  spectrum 
which  can  be  modelled  by  an  equivalent  circuit,  and  identified  by 
nonlinear  complex  parameter  identification. 

In  the  following,  appropriate  impedance  models  for  PEFC  are 
presented.  Using  these  models,  the  electrochemical  behaviour  of 
a  fuel  cell  stack  is  analysed  by  measuring  impedance  spectra  at 
different  operating  conditions  and  by  identification  of  the  model 
parameters. 

2.  Impedance  models 

2.1.  Contributions  to  the  fuel  cell  impedance 

Impedance  spectra  of  fuel  cell  stacks  are  a  superposition  of  var¬ 
ious  intrinsic  loss  factors  and  transient  processes.  Or  in  the  words 
of  an  electrochemist:  “A  multitude  of  fundamental  microscopic 
processes  take  place  throughout  the  cell  when  it  is  electrically 
stimulated  and,  in  concert,  lead  to  the  overall  electrical  response” 
[5].  These  processes  include  the  transport  of  electrons  through 
the  conductors,  the  transport  of  protons  through  the  electrolyte, 
the  transfer  of  electrons  or  ions  at  the  electrode-electrolyte  inter¬ 
faces,  and  kinetic  losses.  In  measured  impedance  spectra  these 
processes  are  not  imaged  separately.  What  can  be  observed  mea¬ 
suring  an  impedance  spectrum  are  dissipative,  capacitive,  and 
inductive  effects.  Finally,  by  interpreting  an  impedance  spectrum, 
these  observations  can  be  assigned  to  dominant  physical  processes. 
Dissipative  effects  can  be  assigned  to  the  ohmic  losses  of  the  elec¬ 
tron  and  proton  conductance.  Capacitive  effects  occur  due  to  the 
double  layers  of  the  electrode-electrolyte  interfaces.  The  dynamics 
of  mass  transport  is  imaged  in  an  impedance  spectrum  as  a  con¬ 
ductive  loop.  Inductive  effects  appear  due  to  the  inductance  of  the 
cables. 

A  difficulty  in  interpreting  impedance  spectra  is  that  besides  sin¬ 
gle  physical  effects  the  impedance  spectra  of  half-cells  superpose 
to  the  impedance  spectra  of  a  single  cell  and  the  impedance  spec¬ 
tra  of  single  cells  to  the  impedance  spectrum  of  a  fuel  cell  stack. 
Reversely,  impedance  spectra  of  fuel  cell  stacks  always  incorporate 
valuable  information  on  the  homogeneity  of  the  single  cells. 

The  purpose  of  an  impedance  characterisation  of  a  fuel  cell 
stack  cannot  be  to  reproduce  all  microscopic  effects.  The  purpose 
of  the  following  impedance  analysis  is  to  determine  the  essential 


properties  and  dominant  loss  factors  of  the  fuel  cell  stack,  their 
interrelations  and  their  dependence  on  controllable  variables  as 
humidification,  conversion,  and  current  density. 

2.2.  Impedance  modelling 

Impedance  models  consist  of  a  mathematical  structure  and 
therein  incorporated  parameters  which  are  constant  for  a  con¬ 
stant  operating  point  but  can  vary  for  different  operating  points. 
Impedance  models  can  be  classified  as  process  models  and 
measurement  models  [6].  Process  models  are  based  on  the 
mathematical  description  of  physical  phenomena  and  therefore 
incorporate  physically  meaningful  parameters  which  can  be  identi¬ 
fied  by  measurement  and  optimisation.  In  contrast,  measurement 
models  are  mathematical  structures  with  shapes  that  fit  to  mea¬ 
sured  impedance  spectra.  Their  parameters  do  not  directly  have  a 
physical  meaning.  But,  by  identifying  the  parameters  for  different 
operation  conditions  physical  phenomena  can  be  assigned. 

The  goal  of  the  following  impedance  modelling  is  to  find  a  model 
with  a  concise  mathematical  structure  which  regards  the  dominant 
effects  and  reflects  the  essential  properties  of  the  electrochemical 
system.  To  enable  an  interpretation  of  the  model  parameters  and 
to  guarantee  a  comparability  of  identified  impedance  spectra,  the 
impedance  model  should  include  a  minimum  number  of  parame¬ 
ters  and  should  yield  a  minimum  error  in  between  the  modelled 
and  the  measured  spectra. 

Fig.  2  shows  a  typical  impedance  spectrum  of  a  polymer  elec¬ 
trolyte  fuel  cell.  It  is  the  Nyquist  plot  of  a  five  cell  stack  with  an 
active  surface  of  100  cm2,  a  commercially  available  membrane  elec¬ 
trode  assembly  with  Pt/Ru  and  the  gas  diffusion  layer  SGL 10BB.  The 
spectrum  was  measured  at  a  current  density  of  0.5  A  cm-2  with  a 
minimum  frequency  of  /min  =  0.2  Hz,  and  a  maximum  frequency  of 

/max  =  10  kHz. 

The  impedance  spectrum  corresponds  to  the  frequency 
response  of  the  fuel  cell  at  an  operating  point  of  the  IV- 
characteristics.  The  task  of  modelling  the  impedance  is  to  find  a 
mathematical  model  which  has  a  frequency  response  with  the  same 
shape  and  size  or  respectively,  to  find  an  equivalent  circuit  which 
shows  the  identical  electrical  behaviour. 

2.3.  Impedance  elements 

The  basic  elements  for  describing  impedances  are  ohmic  resis¬ 
tances  R,  capacitances  C  and  inductances  L  with  their  complex 
resistances:  Zr  =  R,  Zc  =  (jTuC)_1  and  Zi  =  ']coL.  Consequently,  inter¬ 
cepts  of  the  positive  real  axis  refer  to  ohmic  resistances.  A  diverging 


Fig.  2.  Impedance  spectrum  of  a  PE  fuel  cell  stack,  five  cells,  A  =  100  cm2  at  A,c  =  5 
and  A,a  =  1.43,  ipc  =  0.5  A  cm-2,  rstack  =  550C,  TDP  =  50  °C,/min  =  0.2  Hz,  solid  symbols 
denote  decades  of  frequency. 
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impedance 

lim  ilmiz(jco)\\=+oc  (1) 

CD-?  OG  1  1  J  J 

refers  to  an  inductance  L. 

Impedance  arcs  can  be  described  by  a  combination  of  basic  ele¬ 
ments  or  by  distributed  elements.  The  following  two  combined 
elements  are  typical  measurement  models. 

2.3 A.  RC-circuit 

An  RC-circuit  is  a  parallel  connection  of  an  ohmic  resistance  and 
a  capacitance.  The  impedance  of  an  RC-circuit 

(2) 

describes  a  semi-circle  in  the  lower  complex  half-plane  which  is 
called  a  capacitive  loop.  As  can  be  seen  in  Fig.  3(a),  the  impedance 
of  the  RC-circuit  approaches  the  real  axes  for  both  limits,  co  ->  0  and 
co^oc,  with  an  angle  to  the  perpendicular  of  0°.  The  product  r  =  RC 
is  the  time  constant  of  the  RC-circuit. 


2.3.2.  RLC-circuit 

The  RLC-circuit  is  the  parallel  extension  of  an  RC-circuit  by  a 
series  connection  of  a  resistor  Rb  and  an  inductance  L.  Its  impedance 


ZrlcQm)  = 


RaRb  +  J  0)RaL 

Ra  +  Rb  +  j coLRaRbC  +  (j \co)  RaL 


(3) 


corresponds  to  a  second  order  frequency  response.  The  extension 
adds  an  inductive  loop  to  the  capacitive  loop  of  the  RC-circuit  as 
shown  in  Fig.  3(b)  [7]. 

The  following  distributed  elements,  the  Nernst  impedance  Zn 
and  the  generalised  Warburg  impedance  Zw,  are  process  models 
describing  diffusion  processes. 


2.3.3.  Nernst  impedance 
The  Nernst  impedance 


ZnO)  =  rN 


tanh(SNV/jft)/D) 

Vi ®~d 


(4) 


with  the  diffusion  coefficient  D  and  the  diffusion  length  8 n 
describes  the  transport  behaviour  of  the  Nernst  diffusion  [8].  The 
static  resistance  of  the  of  the  Nernst  impedance  is 

Rn  =  ZN(0)  =  lim  {ZN(  j©)}  =  ^  (5) 

CD^O  U 


2.3.4.  Warburg  impedance 

The  generalised  Warburg  impedance 


Zw(j^) 


tanh(Rwa^j&>) 


(6) 


with  the  static  resistance  Rw  can  be  regarded  as  a  bounded  constant 
phase  element  [9],  or  as  a  Warburg  diffusion  element  [10],  where 
a  is  a  function  of  the  diffusion  coefficient  and  the  diffusion  length. 
The  Warburg  impedance  approaches  the  origin  for  co  oo  with  an 
angle  to  the  real  axes  of  45°  as  depicted  in  Fig.  3(c). 


2.4.  Impedance  models  of  a  PE  fuel  cell 


In  Fig.  4  different  impedance  elements  are  assigned  to  the 
impedance  spectrum  of  Fig.  2. 

The  intersection  of  the  impedance  spectrum  with  the  real  axis 
for  co  ->  oo  in  Fig.  4  approximately  corresponds  to  the  ohmic  resis¬ 
tance  Rq,  which  is  mainly  due  to  the  proton  transport  resistance 
of  the  membrane.  The  diverging  impedance  for  co  ->  oo  results  from 


Fig.  3.  Impedance  spectra  of  elements,  (a)  Impedance  spectrum  of  an  RC  circuit, 
(b)  Impedance  spectrum  of  an  RLC  circuit,  (c)  Impedance  spectrum  of  a  Warburg 
impedance. 
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Fig.  4.  Assignment  of  impedance  elements  to  the  impedance  spectrum  of  Fig.  2. 
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the  inductance  of  the  cables,  which  is  not  an  integral  part  of  the  fuel 
cell.  For  intermediate  frequencies  the  impedance  spectrum  consists 
of  two  arcs,  Z^(jco)  and  Z2(jco),  which  have  to  be  examined  in  more 
detail  later.  In  general,  the  impedance  of  a  PEFC  can  be  depicted  as 
the  series  connection  of  the  ohmic  resistance  R two  impedances 
and  Z2  for  the  two  arcs  in  the  Nyquist  plot  and  an  inductance  L. 
The  resulting  fuel  cell  impedance  is 

ZfcCj®)  =  %  +  Z\  (j  CO)  +  Z2(jco)  +  ZL(j(o)  (7) 

The  intersection  of  the  impedance  spectrum  with  the  real  axis 
for  oo  ->  0  corresponds  to  the  static  resistance 


Rz  =  Zpc(O)  =  lim{ZFC(j&>)}  (8) 

o 

which  is  the  limit  value  of  the  fuel  cell  impedance  ZFC(jco)  and  can 
be  calculated  as  the  sum  =  Rq  +  Ri  +  R2,  when  R i  and  R2  are  the 
limit  values  of  Z\  and  Z2  for  co  ->  0. 

In  the  following  four  impedance  models  of  the  PEFC  are  pre¬ 
sented.  These  models  have  the  ohmic  resistance  and  the  cable 
inductance  in  common  and  can  be  distinguished  by  the  ele¬ 
ments  describing  the  two  arcs  in  the  Nyquist  plot.  Hence,  the 
nomenclature  depends  on  the  impedance  elements  describing  the 
impedance  arcs. 

The  structure  of  the  impedance  models  can  be  used  for  the 
whole  operating  range  of  a  fuel  cell.  Identified  parameters  are  only 
valid  in  a  certain  neighbourhood  of  an  operating  point. 

The  first  two  models,  RC-RC  and  Randles,  are  well  known  in 
literature.  The  models  RC-W  and  RLC-W  were  designed  to  fit  the 
impedance  spectra  of  the  examined  fuel  cell. 

2.4.1.  RC-RC  model 

In  the  RC-RC  model  depicted  in  Fig.  5(a),  the  losses  of  the  elec¬ 
tron  and  proton  conductance  are  described  by  an  ohmic  resistance 
and  the  two  loops  of  the  impedance  spectrum  each  by  an  RC-circuit 
[11,12].  The  impedance  of  the  model  is 

ZrcrcCM  =  Rq  +  1  +  j  col  (9) 

1  +)(DT\  1  +J  00X2 

with  the  time  constants  ti=RiCi  and  r2=R2C2  of  the  two 
RC-circuits.  The  series  connection  of  an  ohmic  resistance  and  RC- 
circuits  is  also  known  as  Voigt  electrical  circuit  [6].  The  parameter 
vector 

0rcrc  =  [fli?  R\  Ri  Ct  C2  l]T  (10) 

comprises  the  six  parameters  of  the  model.  The  static  resistance 
equals  the  sum  of  the  resistances. 

This  two-time-constants  model  has  its  limitations,  e.g.,  when 
one  limiting  process  is  diffusion  based.  Depressed  semi-circles  and 
angles  to  the  perpendicular  to  the  real  axes  different  from  0°  cannot 
be  described  correctly.  The  RC-RC  model  can  be  regarded  as  the  first 
approximation  of  the  impedance  of  a  PEFC. 

2.4.2.  Randles  model 

The  Randles  equivalent  circuit  in  Fig.  5(b)  describes  the  response 
of  a  single-step  charge-transfer  process  with  diffusion  of  reactants 
to  the  interface  [10]  and  is  often  found  in  literature  [3,9].  The  equiv¬ 
alent  circuit  consists  of  an  ohmic  resistance  R^,  a  double  layer 
capacity  Cdl  parallel  to  a  series  connection  of  a  charge  transfer  resis¬ 
tance  Rct  and  a  diffusion  impedance  ZN.  With  the  cable  inductance 
the  impedance  is 

^Randles(j^)  =  Rq  +  [j^Ql  +  (^ct  +Zn(J&0)  ^  +  jcoL  (11) 

The  parameter  vector 

^Randles  =  [Rfl  Rct  rN  5n  D  Cdl  I]T  (12) 


Cl  c2 

(a)  Model  RC-RC 


c 


(b)  Model  Randle 


Rn  r 

— dh 


Ri 


Zw  L 


(c)  Model  RC-W 


Rla 


Cl 


Zw  L2 
^ — 0 


(d)  Model  RLC-W 


Fig.  5.  Impedance  models  of  a  PEFC  in  the  form  of  equivalent  circuits. 


comprises  seven  parameters.  The  static  resistance  is  with  (5) 

^Randles  =  Rq  +  Rct  +  (13) 

The  Randles  model  has  the  ability  to  describe  depressed  semi¬ 
circles.  As  the  RC-RC  model  it  has  angles  to  the  perpendicular  to 
the  real  axes  of  0°. 

2.4.3.  RC-W  model 

The  RC-W  model  in  Fig.  5(c)  consists  of  an  RC-circuit  for  the 
first  impedance  arc  in  the  Nyquist  plot  for  low  frequencies  and  a 
Warburg  impedance  for  the  second.  The  impedance  is 

D 

Zrcw(j’&0  =  Rq+  — 7—  +  Zw(j&0  +  )(oL  (14) 

1  +](DT  1 

The  parameter  vector  of  the  RC-W  model 

^RCW  —  [Rn  Rl  Rw  Cl  CL  L]  (15) 

consists  of  six  parameters.  As  at  the  RC-RC  model  the  static  resis¬ 
tance  equals  the  sum  of  the  resistances. 

The  RC-W  model  is  able,  due  to  the  Warburg  impedance,  to 
describe  depressed  semi-circles  and  has  an  angle  to  the  real  axes 
for  co  ^  00  of  45°. 

2.4.4.  RLC-W  model 

If  inductive  effects  occur  in  the  impedance  spectrum  for  small 
frequencies,  the  impedance  element  describing  the  first  arc  in  the 
impedance  spectra  can  be  extended  by  an  inductance.  Hence  the 
RLC-W  model  in  Fig.  5(d)  consists  of  an  RLC-circuit  for  the  first  arc 
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and  a  Warburg  impedance  for  the  second: 

ZrlcwCj’&O  =  Rq  +  ZRLc(jw)  +  Zw(jtu)  +  j  coL  (16) 

The  RLC-W  model  comprises  eight  parameters 

^RLCW  —  [Rn  Rla  Rib  ^W  Q  h  Oi  L2  ]  (17) 

The  static  resistance  of  the  RLC-W  model 

KrLCW  =  Rn  +  n  la,  nb  +  (18) 

Rla  +R\b 

is  the  sum  of  the  ohmic  resistance,  the  Warburg  resistance  and  the 
parallel  connection  of  the  resistances  Ra  and  Rb. 

3.  Electrochemical  impedance  spectroscopy 

The  electrochemical  impedance  spectroscopy  is  an  established 
and  recognised  frequency  domain  method  for  the  characterisation 
of  electrochemical  systems.  It  is  extensively  described  in  litera¬ 
ture  [12-14]  and  widely  used  for  the  characterisation  of  fuel  cell 
impedances  [3,4,15-18]. 

The  procedure  of  the  EIS  is 

1.  harmonic  exctiation  of  the  system  at  defined  frequencies, 

2.  mesurement  of  the  system  response, 

3.  determination  of  the  impedance  at  the  current  frequency, 

4.  parameter  optimisation  of  the  impedance  model. 

Modelling  and  parameter  identification  are  integral  parts  of  the 
method.  The  steps  1-3  are  conducted  for  every  single  frequency. 

In  the  application  of  the  method  for  the  analysis  of  fuel  cells,  the 
fuel  cell  current  at  a  constant  operating  point  is  superimposed  by 
a  harmonic  oscillation 

he  —  Jop  +  ai  sin  (cot)  (19) 

with  the  constant  current  /op,  the  amplitude  a*  and  the  angu¬ 
lar  frequency  co.  The  current  acts  as  system  excitation;  the 
resulting  fuel  cell  voltage  VFC  as  system  response.  With  its  non¬ 
linear  IV-characteristics  a  fuel  cell  generally  does  not  fulfil  the 
Kramer-Kronig  conditions  which  are  linearity,  causality,  stability, 
and  finiteness  [9].  To  meet  the  linearity  conditions  the  ampli¬ 
tude  az  of  the  excitation  signal  needs  to  be  small.  On  the  other 
hand,  the  resulting  voltage  signal  needs  a  minimum  signal  to 
noise  ratio.  Here,  the  amplitude  of  the  current  oscillation  is  cho¬ 
sen  to  0.01  A  cm-2.  Consequently  to  the  nonlinearity  of  the  system, 
impedance  spectra  of  fuel  cells  are  only  valid  in  a  neighbourhood 
of  the  operating  point. 

The  advantage  of  the  electrochemical  impedance  spectroscopy 
is  the  high  accuracy  of  identified  impedance  models  and  its  relia¬ 
bility.  The  main  shortcomings  are  the  long  measurement  time  and 
the  expensive  measuring  setup. 

3.1.  Hardware  setup  and  execution  of  the  test 

For  the  measurement  of  the  impedance  spectra  the  Electro¬ 
chemical  Workstation  IM6  and  the  electronic  load  EL300  from 
ZAHNER-Elektrik  are  used.  With  water  cooling  the  electronic  load 
has  a  maximum  power  of  300  W,  a  current  range  of  0-100  A  and  a 
frequency  range  of  10  mHz  to  10  kHz.  For  the  following  analyses  the 
frequency  is  logarithmically  increased  from  0.2  Hz  to  10  kHz  with 
10  steps  per  decade. 

Prior  to  an  impedance  measurement,  the  controllable  states  are 
set  to  their  operating  points  and  it  is  waited  till  the  voltage  sig¬ 
nal  reaches  a  stationary  level  and  stays  there  for  several  minutes. 
The  waiting  time  is  5-15  min,  depending  on  the  variations  of  the 
operating  points. 


3.2.  Parameter  identification 

The  task  of  the  optimisation  algorithm  is  to  find  the  optimal  set 
of  parameters  which  minimises  the  error  between  the  measured 
impedance  dataZFC  and  the  model  output  ZFC  mod(0).  The  according 
objective  function  of  the  parameter  optimisation  is  defined  as 


k  k 

(20) 

with  the  errors  of  the  real  and  the  imaginary  part  for  every  fre¬ 
quency  cok 

eR,k  =  Re  {ZFC  (j^k)}  -  Re  {ZFC, mod  (j &>k,  0)  } 

(21) 

and 

e\,k  =  Irn  {ZFC  (j cok)}  -  Im  {ZFQmod  (j cok,  0) } 

(22) 

and  the  corresponding  weighting  factors  gR  =  1  and  g\  =  2. 

As  optimisation  algorithm  a  coupled  evolutionary  algorithm 
[19-21]  with  the  deterministic  Nelder-Mead  method  [22]  is  used. 
The  hybridisation  of  the  stochastic  with  the  deterministic  method 
provides  the  possibility  to  combine  the  advantages  of  both  types 
of  methods  [23].  The  evolutionary  algorithm  performs  the  task  of 
the  global  search.  The  Nelder-Mead  method  accelerates  the  con¬ 
vergence  towards  a  minimum.  Ultimately  the  convergence  to  the 
global  minimum  cannot  be  proven,  but  repeated  experiments  have 
yielded  reliable  and  reproducible  results. 

4.  Analysis  of  PE  fuel  cell  impedance  spectra 

The  model-based  analysis  of  the  electrochemical  behaviour  is 
conducted  using  a  fuel  cell  stack  with  five  cells,  an  active  surface  of 
100  cm2,  a  commercially  available  membrane  electrode  assembly 
with  Pt/Ru  and  the  gas  diffusion  layer  SGL 10BB.  The  fuel  cell  is  fed 
with  hydrogen  and  compressed  air,  which  is  humidified  by  a  dew 
point  humidifier.  All  impedance  spectra  are  measured  at  ambient 
pressure. 

4.1.  Comparison  of  the  impedance  models 

4.1.1.  Testl 

At  a  constant  operating  point  the  impedance  spectrum  of  the 
fuel  cell  stack  is  measured.  On  the  basis  of  the  measured  impedance 
spectrum  the  parameters  of  the  models  discussed  in  Section  3  are 
identified  and  compared.  The  operating  point  is  determined  by  a 
current  density  of  iFc  =  0.5  A  cm-2,  an  excess  ratio  of  oxygen  and 
hydrogen  of  Aq2  =  3.33  and  Ah2  =  1.25,  a  fuel  cell  temperature 
of  Tfc  =  55°C  and  a  temperature  of  the  dew  point  humidifier  of 
TDp  =  55  °C. 

4.1.2.  Observation 

Fig.  6(a)  and  (b)  shows  the  Nyquist  plot  of  the  measured 
and  modelled  impedance  spectra.  Fig  7(a)  and  Fig.  7(b)  shows 
the  squared  error,  ek  =  erk  +  eik,  of  the  modelled  spectra  at  every 
frequency.  All  four  models  are  able  to  describe  the  measured 
impedance  spectrum. 

Both  models  with  the  Warburg  impedance,  RC-W  and  RLC-W, 
yield  the  lowest  values  of  the  objective  function  (20).  The  diagrams 
of  the  squared  errors  show  that  the  models  RC-RC,  Randles  and 
RC-W  yield  similar  errors  for  small  frequencies.  Solely  the  RLC-W 
is  able,  due  to  the  additional  inductance,  to  reduce  the  error  at 
small  frequencies.  At  high  frequencies  the  models  with  the  Warburg 
impedance  show  lower  errors,  because  the  Warburg  impedance 
leads  to  an  angle  to  the  perpendicular  on  the  real  axis  of  45°. 
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Fig.  6.  Test  1:  comparison  of  fuel  cell  impedance  models,  (a)  Models  RC-RC  and 
Randles,  (b)  Models  RC-W  and  RLC-W. 


Table  1 

Comparison  of  the  identified  resistances  and  the  values  of  the  objective  function  of 
the  models  RC-RC,  Randles,  RC-W  und  RLC-W.  Here,  R i  and  R2  are  the  limit  values 
of  the  first  and  the  second  impedance  arc. 


Model 

Rq  (mft) 

Ri  (m£2) 

R2  (m£2) 

Rs  (mfi) 

J 

RC-RC 

4.76 

12.9 

7.14 

24.8 

6.84  e-6 

Randles 

4.72 

12.3 

7.74 

24.8 

5.66  e-6 

RC-W 

4.04 

12.3 

8.52 

24.8 

2.82  e-6 

RLC-W 

4.29 

4.82 

15.2 

24.3 

1.48  e-6 

The  static  resistance  is  identically  identified  for  the  models 
RC-RC,  Randles  and  RC-W  (see  Table  1).  The  RLC-W  yields  a  little 
smaller  static  resistance.  The  single  resistances  vary.  Especially  the 
ohmic  resistance  Rq  is  overestimated  by  the  models  RC-RC  and 
Randles  (Tables  2  and  3) 

4.1.3.  Results 

The  RC-W  model  meets  best  the  requirements  for  the 
impedance  modelling.  It  is  able  to  fit  the  measured  impedance  spec¬ 
trum,  yields  a  very  small  error  over  the  whole  frequency  range,  and 
incorporates  a  minimum  number  of  parameters.  Therefore,  it  will 
be  used  for  the  following  analysis  at  different  operating  points. 

4.2.  Impedance  spectra  along  the  IV- characteristics 
4.2.1.  Test  2 

At  defined  operating  points  along  the  IV-characteristics  the 
impedance  spectra  of  the  fuel  cell  stack  are  measured  and  identified 
using  the  RC-W  model.  The  controllable  states  are  set  to  A.q2  =3.33, 
Ah2  =  1.25,  TpC  =  55  °C  and  TDP  =  55°C. 


X  10 7 


Fig.  7.  Test  1 :  error  structure  of  the  comparison  of  the  fuel  cell  impedance  models 
(a)  models  RC-RC  and  Randles,  (b)  Models  RC-W  and  RLC-W. 


4.2.2.  Observation 

Fig.  8(a)  and  (b)  shows  the  impedance  spectra  of  the  identified 
models  at  the  defined  operating  points  along  the  IV-characteristics. 
The  ohmic  resistance  which  almost  coincides  with  the  intercept  of 
the  impedance  spectra  with  the  real  axis  for  high  frequencies  is 
equal  for  all  operating  points.  The  over-all  resistance  decreases 
from  10  to  40  A,  and  from  there  on  increases  till  the  maximum  cur¬ 
rent  of  80  A.  As  can  be  seen  in  Fig.  8(a),  for  small  currents  the  second 
impedance  arc  at  high  frequencies  dominates.  In  contrast,  Fig.  8(b) 
shows  that  for  high  currents  the  first  impedance  arc  dominates. 

Fig.  9  shows  a  diagram  of  the  single  parameters  in  dependency 
on  the  fuel  cell  current.  It  emphasises  the  observation  that  the 
ohmic  resistance  stays  constant  at  constant  humidification  for  all 
currents.  The  resistance  Ri  of  the  RC-circuit  which  describes  the 
first  impedance  arc  for  low  frequencies  increases  with  increas¬ 
ing  current  density.  Whereas  the  resistance  Rw  of  the  Warburg 
impedance  decreases. 


Table  2 

Identified  resistance  values  of  the  RC-W  model  at  varied  hydrogen  excess  ratios. 


Table  3 

Identified  resistance  values  of  the  RC-W  model  at  varied  oxygen  excess  ratios. 
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40A - 60A  - 70A  .  80A 


Fig.  8.  Test  2:  Impedance  spectra  at  operating  points  along  the  IV-characteristics. 
(a)  1=  10-30  A.  (b)  /=40-80  A. 


4.2.3.  Results 

The  structure  of  the  RC-W  model  as  measurement  model  is  able 
to  describe  the  impedance  spectra  over  the  whole  current  den¬ 
sity  range.  Thereby,  the  identified  parameters  show  a  functional 
dependency  on  the  fuel  cell  current.  The  decreasing  resistance  of 
the  Warburg  impedance  and  the  increasing  resistance  of  the  RC- 
circuit  with  increasing  current  density  link  the  parameters  to  the 
overvoltages  of  the  static  IV-characteristics.  Thus,  the  RC-circuit  is 
dominated  by  transport  properties  and  the  Warburg  impedance  by 
charge  transfer  phenomena.  This  result  is  in  contradiction  to  the 
physical  origin  of  the  Warburg  impedance;  but  it  is  in  accordance 
to  the  results  and  assignments  of  the  impedance  arcs  in  [24].  Hence, 
the  mathematical  structure  of  the  Warburg  impedance  is  used  as  a 
measurement  model  to  describe  a  different  physical  phenomenon 
than  the  diffusion. 


10  20  30  40  50  60  70  80 

/(A)  -> 

Fig.  9.  Test  2:  Identified  resistances  along  the  IV-characteristics. 


(a)  -  1,33  -  1,25  -  1,18  .  1,11 


(b)  -  4,0  -  3,33  -  2,86  .  2,5 


Fig.  10.  Test  3 :  Impedance  spectra  at  a  varied  excess  ratio,  (a)  Hydrogen,  (b)  Oxygen. 

The  assignment  of  the  impedance  arcs  shall  be  further  analysed 
in  the  following  tests. 

4.3.  Impedance  spectra  at  a  varied  conversion  of  oxygen  and 
hydrogen 

4.3.1.  Test  3 

In  the  first  test  series  the  excess  ratio  of  hydrogen,  which  is 
defined  as  the  inverse  of  the  conversion,  is  varied  between  Ah2  = 
1.33  and  A,h2  =  1.11  whilst  the  excess  ratio  of  oxygen  is  kept  con¬ 
stant  at  Xq2  =  3.33.  In  the  second  test  series  the  excess  ratio  of 
oxygen  is  varied  between  Xq2  =  4  and  Xq2  =2.5  whilst  the  excess 
ratio  of  oxygen  is  kept  constant  at  A.H2  =  1.25.  The  controllable 
states  are  set  to  iFC  =  0.5  A  cm-2,  TFc  =  55  °C  and  TDP  =  55  °C. 

4.3.2.  Observation 

Fig.  10(a)  shows  the  impedance  spectra  at  a  varied  hydrogen 
excess  ratio.  The  intersection  with  the  real  axis  and  therewith 
the  ohmic  resistance  as  well  as  the  second  impedance  arc  for 
high  frequencies  stay  unchanged  for  all  different  mass  flows.  The 
first  impedance  arc  increases  slowly  with  decreasing  excess  ratio. 
Fig.  10(b)  shows  impedance  spectra  at  a  varied  oxygen  excess  ratio. 
Again,  the  second  arc  stays  constant.  The  first  arc  increases  rapidly 
with  decreasing  excess  ratio. 

4.3.3.  Results 

The  increasing  first  impedance  arc  with  decreasing  excess  ratio 
is  in  accordance  to  the  results  of  the  preceding  test  series,  since  a 
variation  of  the  excess  ratio  primarily  has  an  effect  on  the  transport 
properties.  A  reduced  incoming  mass  flow  of  the  educts  leads  to  a 
lower  partial  pressure  at  the  active  surface,  to  a  reduced  concen¬ 
tration  gradient  and  therewith  to  higher  transport  losses.  A  varied 
excess  ratio  has  less  effect  at  the  anode  side  with  hydrogen  than 
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Fig.  11.  Test  4:  Impedance  spectra  at  a  varied  humidification. 


at  the  cathode  with  the  oxygen  in  the  air.  At  Ah2  >  1.11  they  are 
approximately  negligible.  The  constant  ohmic  resistance  implies 
that  the  fuel  cell  is  sufficiently  humidified  for  all  excess  ratios. 

4.4.  Impedance  spectra  at  a  varied  humidification 

4.4.1.  Test  4 

The  dew  point  temperature  in  the  humidifier  is  stepwise 
increased.  Thereby,  the  relative  humidity  of  the  supplied  air 
increases.  The  controllable  states  are  set  to  iFC  =  0.5Acm_  2,  Aq2  = 
3.33,  A.h2  =  1 .25,  and  TFC  =  55  °C. 

4.4.2.  Observation 

Fig.  11  shows  the  impedance  spectra  at  a  variation  of  the  dew 
point  temperature  between  TDP  =  28  °C  and  TDP  =  55  °C.  The  last 
value  is  equal  to  the  constant  fuel  cell  temperature.  The  identi¬ 
fied  impedance  spectra  marginally  change  their  shape  and  size.  By 
varying  the  humidification  they  are  solely  horizontally  shifted  due 
to  a  decreased  ohmic  resistance  with  an  increased  water  content 
of  the  membrane.  This  result  is  also  reflected  by  the  values  of  the 
parameters  in  Table  4.  With  increasing  humidity  in  the  fuel  cell 
the  transport  hindrances  and  therewith  the  resistance  of  the  RC- 
circuit  increase,  especially  when  the  dew  point  temperature  of  the 
incoming  air  is  higher  than  the  fuel  cell  temperature. 

4.4.3.  Results 

The  observations  emphasise  the  assumption  that  the  ohmic 
resistance  of  the  impedance  models  is  mainly  due  to  the  proton 
transport  through  the  polymer  membrane  with  its  water  content 
dependent  conductivity.  The  increasing  resistance  of  the  RC-circuit 
with  increasing  humidity  again  supports  the  result  that  the  first 
impedance  arc  at  low  frequencies  is  dominated  by  transport  prop¬ 
erties. 

4.5.  Limitations  of  the  impedance  models 

The  RC-W  model  used  in  the  preceding  tests  is  well  qualified  to 
describe  the  dynamic  behaviour  of  the  fuel  cell  stack  at  very  differ¬ 
ent  operating  conditions  in  the  regarded  frequency  range,  fe  [0.2; 
10e3  ]  Hz.  For  frequencies  beyond  the  regarded  range  no  guaranteed 

Table  4 

Identified  resistance  values  of  the  RC-W  model  at  a  varied  humidification. 


statements  can  be  made.  If  an  inductive  loop  appears  for  frequen¬ 
cies  lower  than  0.2  Hz  instead  of  the  RC-W  model  the  extended 
RLC-W  models  has  to  be  used. 

Besides  the  limitation  of  the  model  to  the  regarded  frequency 
range,  the  applicability  is  limited  to  the  measured  current  range, 
JFC  €  [0;  80]  A.  Statements  about  the  electrochemical  behaviour  at 
higher  current  densities  are  not  guaranteed.  But,  the  tests  at  very 
low  stoichiometries  indicate  that  the  model  is  suited  to  describe  the 
electrochemical  behaviour  at  very  high  current  densities,  especially 
the  occurring  transport  limitations. 

5.  Conclusion 

A  simple  impedance  model  is  found  with  a  concise  mathemati¬ 
cal  structure,  a  minimum  number  of  six  parameters  and  a  minimum 
error  to  the  measured  impedance  spectra.  It  is  well  qualified  to 
describe  all  measured  impedance  spectra  at  very  different  oper¬ 
ating  conditions.  By  the  identification  of  the  model  parameters  at 
different  test  series  physical  phenomena  were  assigned.  The  iden¬ 
tified  values  of  the  parameters  at  different  operating  conditions 
provide  the  possibility  to  determine  the  current  state  of  the  fuel 
cell,  to  distinguish  the  major  loss  terms  and  therewith  to  give  an 
insight  into  the  fuel  cell. 

The  impedance  model  with  the  identified  parameters  and  their 
physical  interpretation  are  useful  for  engineering  applications,  e.g., 
to  estimate  the  membrane  resistance  and  for  humidification  con¬ 
trol,  to  observe  the  changes  of  internal  losses  of  fuel  cells  over  time 
without  the  need  for  additional  sensors,  especially  for  lifetime  and 
degradation  testing.  The  model  and  an  appropriate  transformation 
into  time  domain  can  be  used  for  online  diagnostic  methods  [25,26]. 
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